Meron ground states of quantum Hall droplets 
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We argue that topological meron excitations, which are in a strong coupling phase (bound in pairs) 
in infinite quantum Hall ferromagnets, become deconfined in finite size quantum Hall systems. 
Although effectively for larger systems meron energy grows with the size of the system, when 
gyromagnetic ratio is small meron becomes the lowest lying state of a quantum Hall droplet. This 
comes as a consequence of the many-body correlations built in the meron construction that minimize 
the interaction energy. We demonstrate this by using mean field ansatzes for meron wave function. 
The ansatzes will enable us to consider much larger system sizes than in the previous work [A. 
Petkovic and M.V. Milovanovic, PRL 98, 066808 (2007)], where fractionalization into merons was 
introduced. 

PACS numbers: 73.43.Cd, 73.21.La 



I. INTRODUCTION 

A quantum dot (QD) in high magnetic fields^ so- 
called quantum Hall droplet^ represents an exciting 
playground for correlation effects in interacting elec- 
tron systems. Some possible effects can be found 
in transport measurements that detect oscillations in 
magnetoconductance 3 - in the interval 2 > v > 1 of filling 
factors. The associated minima of the current amplitude 
are most completely understood taking into account a 
tendency of the system to find itself in depolarized, highly 
correlated ground states at some fractions in between de- 
spite Zeeman cost<^&£ 

In this paper we will introduce closely related depolar- 
ized states as meron topological excitations of the v = 1 
completely polarized droplet, therefore extending topo- 
logical models to a finite system. They will represent 
relevant lowest lying - ground state configurations of the 
droplet when gyromagnetic ratio is small. 

Whenever we think about quasiparticles in small sys- 
tems we are skeptical about their existence or clear cut 
description that we can find in infinite systems. That is 
even more true if we deal with quasiparticles which quan- 
tization can be based on topological considerations like 
skyrmions in quantum Hall (QH) ferromagnets £• Nev- 
ertheless, as we will argue here, topological objects like 
merons, a meron is a half of skyrmion^ can exist as low- 
est lying states of a quantum Hall droplet. Therefore 
fractionalization and quantization, characteristic to QH 
systems, may persist even in small spin unpolarized sys- 
tems bringing topological objects to their description. 

The spin characterization of QD states at realistic Zee- 
man coupling is an open problem* 1 ^ A very thorough un- 
derstanding of the completely polarized case exists that 
includes a description of the QD states by the way of a 
vortex quantization i 11 ' 12 : 13 ' 14 The introduction of higher 
angular momentum states or the increase of magnetic 
field with respect to the most compact, maximum den- 
sity configuration of a dot is followed by vortex appear- 
ance inside the dot. On the other hand there is a need to 



introduce a classification to various (partially polarized) 
states at arbitrary or realistic Zeeman coupling^ 

In the previous work, Ref. [la . by one of the present 
authors, the Coulomb interaction problem of small quan- 
tum Hall droplets (dots) with N = 4 and N = 6, N is 
the number of electrons, in the limit of zero Zeeman cou- 
pling, was studied by exact diagonalization, in the lowest 
Landau level approximation. It was shown that the low- 
est lying states of these small quantum Hall droplets can 
be described and classified as states of merons. 

In this work we address the question of the existence 
of meron ground states in large quantum Hall droplets, 
N ~ 20, that exact diagonalizations can not reach. At 
the same time, by extrapolation, we will be able to es- 
timate the size of the quantum Hall droplet, N ~ 100, 
at which meron confinement takes place i.e. when we 
can not expect meron ground states. We will be able to 
do this by using a model wave function that describes a 
meron of arbitrary winding number positioned at the cen- 
ter of a droplet. In this sense our approach is variational, 
takes into account a mean field description of a meron, 
and compares it energetically first with the ground state 
- maximum density droplet (MDD) spin polarized con- 
figuration, but also with other meron configurations - of 
different winding number. 



II. DESCRIPTION OF MERON GROUND 
STATES 

At the heart of the theory of the quantum Hall systems 
at effective small gyromagnetic ratio, and filling factor 
v = 1, i.e. quantum Hall ferromagnets, is the commen- 
suration of spin and charge deviations from the ground 
state values^ Exchange interaction prefers smooth tum- 
bling of spins that follows changes in the charge distri- 
bution. Topology plays an important role in the the- 
ory of infinite size quantum Hall ferromagnets." After 
the spontaneous symmetry breaking the polarization of 
the ground state is in a definite direction and fixes the 
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boundary condition at infinite radius. Due to then pos- 
sible mapping between real and internal (spin) spaces, 
the skyrmion excitations with nontrivial value of topo- 
logical charge, and the same value of electric charge im- 
plied by the commensuration, can be identified. This 
classical (nonlinear-cr model) skyrmion solution finds a 
concrete quantum-mechanical realization in the follow- 
ing construction 9 



JV 



n : n< 



Zj)- 



(1) 



i<3 



We omitted Gaussian factors and z's are 2D coordinates 
of N electrons in the lowest Landau level (LLL) in units of 
the magnetic length, I b = ^jTic/eB. In this construction 
the coordinates (in the orbital part) participate in the 
spinor part. The complex number A denotes the size of 
the skyrmion signifying the characteristic length when 
the spin, which is around the center pointing up, starts 
gradually pointing down as it does at infinity ( \z\ — > 
oo). The parameter A is finite , determined by the ratio 
between the strength of the interaction and the Zeeman 
coupling. 

Now we specialize to the case of a droplet. We sim- 
plify the matters assuming that the gyromagnetic factor 
is zero. Thus we are in a scale invariant situation where 
A may be function only of the size (JV) of the system. 
We take A — * c = \/N jil and consider constructions for 
merons of the following form 



N 



(2) 



i<3 



where n is a positive integer. The mean angular momen- 
tum M can be easily (numerically) calculated and with 
respect to the ground state value, M a = JV(JV — l)/2, we 
find that 



AM = M - M Q =n— 
2 



JV n(n - 1) 



(3) 



approximate very well the calculated AM, the better, 
the larger JV is (Fig. 1). Therefore, for n = 1, the excita- 
tion, in the mean, is the excitation of one half of the flux 
quantum, AM — N being the excitation of (Laughlin's 
quasi) hole of one flux quantum or a vortex. A skyrmion 
is a generalization of a Laughlin quasihole in the case of 
quantum Hall ferromagnets. We may then identify the 
construction with a classical one for a single meron being 
" half of the skyrmion" , and carrying flux of one half of 
one flux quantum (hc/e). Also its magnetization, from 
pointing up in the center, gradually transforms into one 
in the plane (S z = 0) on the edges resembling the half 
of the skyrmion^ In the case of a skyrmion the magne- 
tization would proceed to transform into configuration 
that points down on the edges with equal magnitude but 
opposite direction with respect to the configuration at 
the center. For a detailed description of the single meron 
(Eq. © with n — 1) see below and Fig. 5. 



AM / JV= n/2 + n (n-1) I IN 
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FIG. 1: The mean angular momentum with respect to the 
ground state value of the constructions in Eq.((2}. The results 
for various system sizes, i.e. number of electrons (JV), are very 
well fitted with Eq.Q. The fit is better the larger the system 



Now we will show the results of the calculations of the 
total (interaction plus confining) energy with respect to 
a suitable Hamiltonian of the states expressed by Eq.@. 

The interaction part of the Hamiltonian we work 
with can be described as a truncated pseudopotential 
interaction^ with only pseudopotentials Vi for I = 0, 1 
nonzero and positive. Here I is the relative angular mo- 
mentum of a pair of electrons in the LLL. Precisely, 



(4) 



i<j 1=0 



where P[j is the projector of a state in the LLL of the 
pair (ij) to the definite relative angular momentum value 
I. We consider Vq much larger than Vj., making the 
state in Eq.Q with n = a unique ground state (up 
to SU(2) rotation) of the system at filling factor one. In 
addition to Vq ^ we choose Vi ^ because we want 
to break energy degeneracy among topological construc- 
tions implied by the hard core (only Vq ^ 0) models For 
an explanation how to implement the model interaction 
with pseudopotentials see Appendix. That meron excita- 
tions (described in mean field as central constructions by 
Eq.©) are indeed lowest lying states of a small droplet 
was shown in Ref. [l6| by exact diagonalization in the case 
of the Coulomb interaction. 

The interaction energy of the ground state we calcu- 
lated taking Vi = 1 in Hi and the results are plotted in 
Fig. 2. We see the dependence of the interaction energy 
per particle on the size - the number of particles in the 
system. The energy can be very accurately fitted with 
the analytical expression in the Figure as an expansion in 
inverse powers of y/N. ( y/N is the measure of the radius 
of the system.) The same analytical expression can be 
found just considering the sizes up to JV = 20. 

The expectation values (interaction energies per par- 
ticle) of the meron states (Eq.© with n — 1,2,3) with 
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FIG. 2: The ground state interaction energy per electron in 
units of V\. The analytical fit is an expansion in inverse pow- 
ers of yN which is the measure of the radius of the system. 
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FIG. 3: The interaction energy per electron with respect to 
the ground state value for the meron states with n = 1,2, 3. 
Continuous lines depict fits to the calculated values. Fits are 
expanded in powers of 1/yN and the leading term in the 
N — > oo limit is positive. This behavior leads to the absence 
of merons in large systems. 



respect to the ground state value are plotted in Fig. 3 
along with analytical fits. Again we take V± = 1 and 
there are no contributions from the first pseudopotential 
Vo because the states in Eq.@ are zero energy eigen- 
states with respect to that pseudopotential. The analyt- 
ical fits are expansions in inverse powers of VN, which is 
the measure of the radius of the system. The expectation 
values come with overall negative sign with respect to the 
ground state value as each excitation represents an infla- 
tion of the volume of the system and therefore increase 
in the average distance and decrease in the interaction 
among particles. Nevertheless, we can find extrapolating 
the fitted analytical expressions to larger number of par- 
ticles and system sizes that the positive first term (small 
but always present and approximately proportional to a 
constant times n) will overcome all other terms in the 
N — ► oo limit. This leads to the conclusion that in the 
thermodynamic limit, AEj ~ N and positive, the merons 
are not present in the excitation spectrum. Like vortices 
in the XY model they are confined in pairs as skyrmion 




FIG. 4: The total energy gain per electron as a function of n 
for N = 19, and for different values of the confining constant 
g. There is no meron ground state for g > 0.032. Here, 
g C i,i — 2,3,4,5 (c for critical) are calculated values of the 
slope of the confining potential when the total energies for 
creation of merons with winding numbers n = i — 1 and n = i 
are equal. Therefore g C i has the meaning of the slope of the 
confining potential when meron with higher winding number 
i enters the droplet. 

excitations. 

In a strong magnetic field and realistic situation with a 
harmonic confining potential, we can model the confining 
part of the Hamiltonian as 

H c = gAM. (5) 

AM is the angular momentum measured with respect 
to the ground state value, Eq.©, and g is a positive 
constant. Then we may expect (due to the positive third 
terms in the expressions in Fig. 3 with nearly n(n — 1) 
behavior) that in the case of the Hamiltonian, 

H = H I + H C (6) 

the system is prone to the instabilities (spin reconstruc- 
tions) described by the meron excitations in Eq. ^ . Thus 
we plotted the total energy of a droplet with N = 19, 
Fig. 4, for a few g c (c - critical) when the total energy 
of the excitation with fixed n equals the energy of the 
one with n + 1. This is a situation when one ground 
state with fixed n becomes less favorable with decreas- 
ing g and substituted with another one with the meron 
number equal to n + 1. Somehow, as an only exception, 
our mean field approach predicts that the n = 2 state re- 
construction precedes the one with n = 1 ( Fig. 4). This 
is very likely an artifact of our choice of the interaction 
potential ( Eq.Q). The implied ground states are only 
possible when the size of the system is not large that 
the first and the second term in the expansions of the 
interaction energies become comparable in size. Approx- 
imately this happens when N > 100 as can be seen from 
the analytical fits in Fig. 3. 

With a view on the theory of the evolution of a quan- 
tum dot with magnetic field in Ref. 2fj in terms of com- 
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pletely polarized states incorporating (quasi)hole excita- 
tions, and its reasonable agreement with experiment,- 22 
the meron evolution that we advocate would occur at 
much smaller gyromagnetic ratio (or stronger interac- 
tion strength). Further merons of higher meron num- 
ber (n — 2 or 3), "giant merons", are highly unlikely in 
larger droplets, N > 5, as described in Ref. 16. This 
will likely cause that implied N/2 periodicity would be 
slightly modified, as entering merons would accomodate 
also on orbitals away, off the center just in the case of 
vortex excitations in the completely polarized case^ 

To illustrate the bulk single meron configuration that 
the construction in Eq. ([2]) for n = 1 represents we plotted 
in Fig. 5 its spin density S(z), and charge density p(z), as 
functions of the 2D coordinate z = \z\ exp(i<y9) in the case 
of N = 20. The analytic expressions (which are plotted) 
for spin density S(z) are 



S z {z) = (N-\z\ 2 )f N (\z\), 
S x (z) = 2VN\z\cos(<p)f N (\z\), 
S y (z) = 2y/N\z\sin(ip)f N (\z\), 

and for charge density p(z) is 

P (z) = (iv + M 2 )MM), 

where 



MM 



\z\ 2m e~ 



^ 27r2 m+1 m! (m + 1 + f ) 

1—0 V * / 



(7) 
(8) 
(9) 

(10) 
(11) 



Therefore we have S 2 (z) = p 2 (z). From Fig. 5 and these 
expressions we can conclude that for \z\ < V^/V we have a 
description of a single meron that is created in the bulk 
of the system. Close to the boundary S(z) vector lies in 
the 2D plane and rotates for 2tt as expected for a meron. 

Now, we must pose a question which quantum me- 
chanical states - eigenstates of M and S z> correspond to 
the "mean-field" or "classical" constructions expressed 
by Eq.@. We will take that each meron classical con- 
struction corresponds to the quantum mechanical state 
that we find by expanding the spinor part of Eq.© with 
eigenvalue of M corresponding to the expectation value 
in Eq.Q. On the other hand, the states that were pro- 
posed for spin (global not only edge) reconstructions in 
Ref. SH are of the following form (s{) k °\C N >. \C N > 
denotes the filled with TV spin | particles LLL - a Vander- 



monde determinant, and S{ = J2m v / w + 1 c, 
an exciton operator with c\ n and c m , the electron cre- 
ation and annihilation operators, and m denoting single- 
particle angular momentum. k a — N/2, if N is even, and 
after a little inspection we can find out that the state de- 
scribed by the formula coincides with our quantum me- 
chanical eigenstate analog of the construction in Eq.© 
with n equal to one. As we now show explicitly, in the 
case of the merons of Eq.@, these constructions carry 
very small S z . The states in Eq.© are eigenstates of 
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FIG. 5: Plotted are for N = 20, in the upper panel, charge 
density p(|z|) with respect to the density in the ground state, 
po(|z|)j in the lower panel, S^dzl) component of the spin den- 
sity along S x {\z\,<p = 0) = y/S$(\z\) + S$(\z\), and in the 
side panel, spin density vector S(|z|) as functions of radius \z\ 
for the meron construction, Eq.([2} with n = 1. A small dip 
in the charge density for \z\ < VN describes the deficit of the 
charge due to the presence of a meron. The amount of the 
missing charge is 0.486e, close to e/2 as expected for a meron. 



AI + nS z — M Q with eigenvalues nN/2H Therefore (S z ) 
in these states is (S z ) = — (n— 1)/2 and n > 1 for merons. 

To find out whether even in small droplets the mean 
field ansatz has a distribution of participating eigenstates 
of M (and S z ) very much peaked around the quantum 
mechanical states with AM given in Eq.Q, we calcu- 
lated these distributions (Fig. 6) in the cases of a droplet 
with N = 17 and N = 20. Although our trial functions 
are not eigenstates of the angular momentum, the calcu- 
lated participation ratios, 



Pi, 



(0 + ion)! • • • (N - 1 + i N -m)\, (12) 



<1 ,...,2JV-1> 



where if = 0, 1 and ^2^=0 ff = ^' are m deed well peaked 
around expectation values, 



AM(N,n) = ^ k =° nkPk . 

V Pi 



(13) 



III. DISCUSSION AND CONCLUSIONS 

The meron physics we described above should be in 
principle detectable in lateral quantum dots, in which 
interaction effects are strong, in their not yet explored 
regime beyond the MDD stated as well just below the 
MDD state where depolarized (S — 0) states were al- 
ready discovered A 
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mentum of N/2, Due to being just vortices of one kind of 
fcrmions, a density of fermions of another kind at their 
core was detected. Because of this phenomenon these 
vortices are called coreless vortices. For a larger droplet 
[N = 20) it was demonstrated using DFT (density func- 
tional theory) and the RWF method that a coreless vor- 
tex enters the droplet (at an angular momentum slightly 
higher than that of MDD). Based on these observations, 
i.e. that coreless vortices as merons show an increase of 
polarization at their core (center), and our studies we can 
conclude that merons are the solutions of the quantum 
Hall droplet in the fully spin rotationally invariant case, 
i.e. in an indistinguishable picture, and correspond to 
these coreless vortices in the distinguishable picture. 



FIG. 6: Participation ratios Pfc/P max of the eigenstates of the 
angular momentum with eigenvalues M — N(N — l)/2 + nfc as 
functions of k = 1, . . . ,N for the meron construction defined 
by Eq.©. 



Besides the study in Ref . KLq of small systems, there are 
studies in Refs. [lQj 24]]25| without the LLL approximation 
that find depolarized states that we can identify as meron 
ground states. They appear at angular momenta, multi- 
plets of N/2, for N even. 

Here we addressed the question of the existence of the 
meron ground states in larger quantum Hall droplets. 
Our conclusion is, as long as the Zeeman term is small, 
the meron ground states are viable solutions of these 
droplets for N smaller than ~ 100. Also we expect the 
appearance of merons in any small enough system of par- 
ticles with the lowest Landau level quantization and a 
degenerate, additional degree of freedom. That may hap- 
pen in graphene structures in the presence of magnetic 
field (with valley degree of freedom) or rapidly rotating 
quantum gases (with spin degree of freedom). 

A previous study, Ref.l26l. of few-electron qnatum dots, 
without the Zeeman term, was based on restricted (condi- 
tional) wave function (RWF) method^ As the method 
underlines explicit fixing of the spin projection of each 
electron, the vortices observed are analogous to those of 
polarized systems. A more detailed and recent study 
in Ref. H?] is based on the same method and discusses 
the problem of two distinguishable species of bosons or 
fermions which rapidly rotate. As usual in these cir- 
cumstances the LLL approximation is applied and the 
fcrmion system is (up to the difference caused by the dis- 
tinguishability) the same as in the case of the quantum 
Hall droplet discussed in this paper. The first vortices 
that enter the MDD come at the increase of angular mo- 
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APPENDIX 

First we consider a fixed spin configuration - compo- 
nent of the total wave function which expectation value 
of the interaction energy we want to calculate: 



,z N ). 



(A.l) 



"J" in general does not have any overall definite symme- 
try property (i.e. it is not antisymmetric or symmetric 
in general). We first have to perform the calculation 
for each component ^ separately and then add contribu- 
tions. We do the calculation by performing the projec- 
tion to angular momentum, 1 = 1, component of wave 
function ^ for each pair of particles Zi and Zj . The com- 
ponent can be extracted by expanding 'J in powers of 
(zi + Zj) and (z« — zj) i.e. 

m,n— 

Cm,n (^1 ) • * * j Zi-t,z i+1 , . . .,Zj-i,Zj+i, . . -,z N ), 

(A.2) 

and taking only m = 1 contribution - wave function in the 
expansion. Because we take Vi = i = 1, the normalized, 
integrated over all coordinates square of the projection 
would be the contribution to the interaction energy from 
the given spin configuration "J and fixed pair of particles. 
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